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Abstract 



I identify a source of ^.QCD/mfiQ corrections to the assumption of quark- 
hadron duality in the application of heavy quark methods to inclusive heavy 
quark decays. These corrections could substantially affect the accuracy of 
such methods in practical applications and in particular compromise their 
utility for the extraction of the Cabibbo-Kobayashi-Maskawa matrix element 
V cb . 



*an abbreviated version of the original JLAB-THY-98-03 entitled "Duality in 
Inclusive Semileptonic Heavy Quark Decay" 
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Although the classic application of heavy quark symmetry is in the exclusive semileptonic 
decays of heavy quarks there has also been substantial work on using heavy quark 
effective theory (HQET) to systematically improve decay predictions for inclusive decays 
of heavy hadrons In these inclusive applications, decays are treated in an operator 

product expansion (OPE) which leads via HQET to a I/uiq expansion in which the leading 
term is free quark decay and I/tuq terms appear to be absent. Although these calculations 
have become very sophisticated f4j||, it is widely appreciated {§-[7j] that there remains a 
basic unproved hypothesis in their derivation: the assumption of quark-hadron duality. It 
is the accuracy of this assumption that I want to call into question here. 

While supposedly of wide validity, recent applications have centered around the hope 
that this approach offers an alternative to the classic exclusive methods for determining V c b, 
and I will accordingly focus most of my remarks on the case b — > civi where both quarks 
are heavy. In inclusive b — > c£i>£ decays, which materialize as B — > X c £ui, about 65% of 
the X c spectrum is known to be due to the very narrow ground states D and D*. The 
relatively narrow sj* = | + states || -D2(2460) and -Di(2420) account for perhaps another 
5% of the rate, and it may be assumed that the remaining rate involves decays to higher 
mass resonances (quarkonia and hydrids) and continua ||. The inclusive calculations predict 
continuous X c spectra which are assumed to be dual to the true hadronic spectrum (see Fig. 
!)• 

A picture like Fig. 1 might lead one to dismiss the duality approximation since the inclu- 
sive spectrum clearly does not meet the usual requirement that it be far above the resonance 
region [|10|1 . I.e., normally the accuracy of quark-hadron duality would be determined by a 



parameter Aqcd/E where the relevent energy scale E is the mean hadronic excitation en- 
ergy Am>x c = mx c — m D- However, as first explained by Shifman and Voloshin [Il],[l2|], this 
is not the expansion variable in this case: duality for heavy-to-heavy semileptonic decays 
sets in at threshold since even as 5m = mi, — m c (and therefore Am^J approaches zero, as 
nib oo the heavy recoiling c quark has an energy much greater than Aqcd so that it is a 
free quark in leading order. In the small velocity (SV) limit, it must therefore hadronize with 



unit probability (up to potential Aqcd/ttiq corrections) as D and D*. This "cannonball" 
approximation is in fact an essential part of the physical basis of the HQET expansion in 
1/rriQ. Thus the question is not whether duality holds in semileptonic heavy quark decays, 
but rather how accurately it holds. 
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Fig. 1: A sketch for b —* c semileptonic decay of the continuous inclusive recoil spectrum 
of the OPE calculations (smooth curve) compared to the known hadronic spectrum (shown 
as individual resonance lines). 



Up to caveats regarding the unknown accuracy of the assumption of duality, the combined 
HQET and OPE methods indicate that inclusive calculations should in fact be accurate up to 
corrections of order h^ CD /m 2 Q. Here I will identify a source of duality- violation which leads 
to Aqcd/i^q corrections for any finite final quark kinetic energy. The problems are revealed 
by considering a Bjorken sum rule which may be viewed as an extension of Shifman- 



Voloshin duality to arbitrary recoils. Bjorken's sum rule guarantees that, as m& — > oo, 
duality will be enforced locally in the semileptonic decay Dalitz plot of rate versus w — 1 and 
Eg (where w = v ■ v' is the usual heavy quark double-velocity variable and Eg is the lepton 
energy). For regions of the Dalitz plot for which w — 1 is not large (and in b — > c decay 
nearly the whole Dalitz plot satisfies this condition), the Bjorken sum rule explicitly relates 



the loss of total rate from the "elastic" sj* = | channels, as the Isgur-Wise function falls, 



to the turn-on of the production of s[ l = ~ and ^ states |14 



d 2 T A 



K(w,E e ) dwdE e 



quark and hadronic total 



elastic 



inelastic 



w-1 



Fig. 2: The exact compensation by inelastic channels of the fall of the elastic rate in the 
linear region as — > oo. 



In particular, in this region the Isgur-Wise function may be taken to be linear: 

1 



£{w) ~ 1 - p\w- 1) = 1 - 



^ + Pdyn 
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and if we define 



J2~ninclusive 
u quark 

dwdEp 
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then [I3l|l4l, as m& — > oo, 
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the Bjorken sum rule guarantees that for fixed r = m c /mb, as m& — > oo inelastic s 



■Kl _ 1" 



and | channels will open up to give a semileptonic rate that exactly and locally compensate 
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in the Dalitz plot the loss of rate from the elastic channels due to p\ yn - This situation is 
sketched in Figure 2; if it were applicable to b — ► c decays, then quark-hadron duality would 
be exact. 

Having established conditions for its validity as m& — > oo, it is easy to see why one should 
be concerned about quark-hadron duality for b — > c decays. For fixed r, w — 1 lies in the 
fixed range from to (1 — r) 2 /2r, and as m& — > oo any given hadronic threshold collapses 
to the point w — 1. However, for finite m& there is a gap in w — 1 in which the rate to the 
elastic | channels falls by Kqcd/^q terms but the potentially compensating excited state 
channels \ + and | + are not yet kinematically allowed. More precisely, if mo** is the mass 
of a generic charmed inelastic state, then = {tub — mo**) 2 would be the threshold in t 
for this state, corresponding to a value of w — 1 in the quark-decay Dalitz plot of 

, A 



J. j.** 



(5) 



2m\ ) m c m c 

where t m = (mg — mp) 2 ~ (mt, — m c ) 2 and A = mo.. — rrio- Since A ~ 500 MeV 
and (w — l) max — 0.6, this region covers more than one third of the Dalitz plot and the 
compensation is very substantially delayed: see Figure 3. Eqs. (5) and (1) show that this 
effect is of order Aqqd/^q, seemingly at odds with the OPE result. 
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Fig. 3: An overlay of the Dalitz plots for B -> D^eP e , B -> D**(2500)ez/ e , and B -> 
D**(3000)ez/ e . The mass is taken as the hyperfine average of the D and D* masses; the 
two D** masses are chosen for illustrative purposes. 



Despite this apparent contradiction, there is actually no inconsistency: the OPE result 
that the leading corrections to the inclusive rate are of order Ag C , D /mg can still be valid as 
derived in the limit of large energy release in the b — > c transition, while Aqcd/^q effects 
can arise for energy releases of the order of Aqcd due to a finite radius of convergence of 
the OPE. The main purpose of this paper is indeed to call attention to this effect. 

The basic issues can be most easily exposed by considering [|J spinless quarks coupled to 
a scalar field <fi of mass //, and by studying the decay b — > afi with weak coupling constant 
g. Differential semileptonic decay rates have a more complex spin structure, but otherwise 
correspond to the case fj, = \ft; total semileptonic rates correspond to a weighted average 
over kinematically allowed fi but, as we shall see below, this averaging does not change the 
essentials of the problem. In our simplified case 

T(b - c0) = |H (6) 

where pji = [(m, — vrif) 2 — /i 2 ] 1//2 [(mj + mj) 2 — fi 2 ] 1 ^ 2 /2rrii is the momentum of <p from the 
two-body decay of mass m; into masses m/ and [i. 

To compare Eq. (6) with a hadronic world (I initially consider a large N c world of narrow 
resonances, but will generalize below), define 

T(B - DM® = ^^(^^)|e (n) (^) B )r (7) 
87rm B v m c mb ' 

where the generalized Isgur-Wise functions depend on 

V D {n) B =P D {n) B /m D { n ) , (8) 

the recoil velocity of the n th excited state of the D meson system, and I have introduced 
some conventional mass factors to explicitly reflect hadronic normalizations. I next introduce 
a "scaled energy release" variable 

m * m b — m c — V> 

T*= A c M , (9) 

where A = vtld** — mr, is the mass gap to the first s[ l = 1~ excited state of the D meson 
system (corresponding to sj' = | + and | + in the physical case where the d has j p = \ ), 
and an order KqcD/ m Q expansion parameter e defined by the expansion 
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\i DB \ 2 = 1 - eT* + 0{t 



(10) 



for small T* {i.e., small charm quark velocities). In the quark model [|T(|,|T7j] one would have 

m d (m b - m c ) 



m h m c 



(11) 



where rrid is the mass of the light spectator antiquark d (or, more generally, of the "brown 
muck"). Defining |£d**b| 2 = J2 me \£,d**b\ 2 (where £™L B is the analog of £db for transitions 
into the lowest s\ l = 1~ excited state with magnetic quantum number mi), we would have 



|^| 2 = e(T*-l) + 0(e 2 



(12) 



from the (spinless) Bjorken sum rule in the limit that it is saturated by the first D** . Since 
in this limit 
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we can obtain a model [HI for 
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3 (T* - l) 3 / 2 

[1 + -e - eT*}6(T*) + e { 6(T* 



(15) 



(16) 



wherein I have shown explicitly the two thresholds at T* = and T* — 1. It is interesting 



to observe that the quark model of Refs. ||lq , |l7|l gives exactly Eq. (16), including the +|e 



term flSfl , as expected [0. I also note that 



1. At T* -> oo, Eq. (16) is of the form 1 + 0(e 2 ) + 0{e/T*) as required by the OPE. 

2. There are no other terms of order 1, e, or eT* possible beyond those shown: a more 
accurate treatment of T(B — > could only generate e 2 , e 2 T*, e 2 T* 2 , ... terms; a more 



accurate treatment of T(B — > D**<p) could only generate e 2 T*, e 2 T* 2 , ... terms; and all 
higher states first make a contribution at order e 2 T* 2 or higher. Conversely, we note that if, 
for example, e 2 T* 2 terms are retained, they must all cancel exactly or the requirements of 
the OPE would be violated as T* — > oo. 

3. As Am = m b -m c ^ 0, flf " -> 1 + 0{ AQCD 2 Am ) as required 0. 

4. Near T^ ax = m b — m c , €T^ ax is in general large. This observation corresponds in the 
usual language of heavy quark symmetry to the statement that the natural scale of the slope 
p 2 of the Isgur-Wise function is of order unity. It is also consistent with the experimental 
observation that |£db| 2 has dropped to less than half its value between zero and maximum 
recoil. Given this, the extension of Eq. (16) to higher orders in T* will require a "conspiracy" 
of the entire spectrum of possible hadronic final states. We may nevertheless use Eq. (16) 
across the full range of T* as an indicator of the Aqcd/ m Q effects arising from the order 1 
and order T* terms in the expansion of R. This corresponds to a "best case" assumption 
that duality is locally perfect for the terms T* n with n > 1. 

Thus, while extreme, this truncation has all the properties required by the OPE and so stands 
as a simple explicit example of the existence of the claimed duality- violating Aqcd/^q 
effects for finite T* . 

It is straightforward to introduce a number of simple variants of this prototypical model. 
The first corresponds to the more realistic case where the Bjorken sum rule is only saturated 
by the full tower of sj* = 1~ resonances so that the second term in Eq. (16) becomes 

Tf^EUT* -t:f 2 9(T* -O (17) 

with J2 n fn — 1 an d t* n being the threshold for channel n. As T* — > oo, these contributions 
automatically cancel the — eT* term from the elastic form factor, and constrain the 0(e) 
correction to give 

R?ll + ... = [1 + \& - tT*]9(T*) + UT* - t* n f/ 2 e(T* - r n ) , (18) 

* J- n 

where 
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Fig. 4: Four resonance models of the approach to duality: (a) i?f ** (with the baseline 
value e = 1/6), (b) i?f" (with e = 1/3), (c) R?+ 2+ .„ (with t* n = n and /„ = (|) n so that 
i* = 2), and (d) R ho . The "baseline" parameters follow from the observations that p 2 ~ 1, 
(w — l) m ax ~ 1/2 and {m b — m c )/A ~ 6. The alternative e = 1/3 corresponds to the case 
Tmax = 3, «.e., to using A e jj = 2 A for the mean location of the sj e = l~ strength. 

As described above, both and -Rf^2+ are still "best case" truncations which assume 
exact cancellations of e 2 T*, e 2 T* 2 , ... terms. While sufficient for the purposes of this study, 
this limitation is easily removed: it is straightforward to recursively "construct duality" to 
create models to the required order in e to any finite order in T*. Consider, for example, 
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([1 + ^ e ]T* 1/2 6(T*) 

^ 1 + r + S e2 + l e3 + ^ 4 + ^^- 1 ) 3/2 ^- 1 ) 
+le 2 [l + 7 -e + ^e 2 + ^e 3 + ...](T* - 2f 2 6{T* - 2) 

+^ 3 [1 + \e + f^e 2 + - 3) 7 ^(T* - 3) 



+ ^e 4 [l + ye + ...](T* - 4) 9 / 2 #(T* - 4) 
+ ^e 5 [l + ...](T* - 5) 11/2 #(T* - 5) + ...) 



where the ellipses denote terms of order t % T* n with 1 < n < 5 and all terms of order t m T* m 
and higher with m > 5. This harmonic-oscillator-like expansion is accurate even at T* = 5 



The three models just introduced are all based on the duality of b — > c0 to a simple 
tower of cJ resonances controlled by the single scale A. Figure 4 shows that the thresholds 
associated with such towers could easily be a source of duality- violating KqcdI m Q correc- 
tions of order 10% in b — > c decays. This must be a cause for concern in comparing inclusive 
calculations with experiment. 

I am even more alarmed by processes which could give a high-mass nonperturbative 
tail to the recoil mass distribution. The hadronization of b — * c<p will not be saturated by 
ordinary quark model cd states even in the large N c limit: hybrid mesons (i.e., states with 
a cd valence structure but with internal gluonic excitation) will also contribute. Such states 
are expected at substantially higher masses than the ordinary quark model states. Moreover, 
their production will not be exhausted until the constituent d antiquark in the D meson has 
been fully resolved into a current quark at high recoil momentum p c >> 1 GeV. For a crude 
estimate of the effects of the delayed onset of these states, I take a simple two-component 
resonance model consisting of "normal" cd resonances with t*^ and cd hybrids with i* hybrid 
substantially larger. If we assume that the latter are responsible for a fraction k of p 2 d , 
then we would have 



up to corrections of order Aq 



QCD 
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Fig. 5: Four examples of the effects of a nonperturbative high mass tail on the recoil mass 
spectrum: (a) R*l y " (with k = 1/10, t*j = 2, and i* hyhrid = 4), (b) R 6 V n (as in (a), but 
with 4^ = 6), (c) i?r (with A = 1/5, r d - = 2, T* in = 2, and s = 4), and (d) iJ™ (as in (c) 
but with s — 6). The values of k and A are based on the model-dependent estimates of Ref. 
0, but are certainly reasonable (e.g., A is a l/N c effect). The illustrative values for i\ ybrid are 
based on the high threshold for hybrids and their presumed "hard" production mechanism. 
The nonresonant spectrum is assumed to have the form p(t*) = ^exp( Tmin s 1 )6(t* — T^ in ) } 
with the choices for s reflecting the assumed persistence of nonresonant contributions to 
invariant masses of order 2 GeV above threshold. 
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I suspect that l/iV c -suppressed nonresonant contributions are an even more serious source 
of delayed compensation of duality. An appropriate model for such continua would be 



normalized spectral function (J^S dt*p(t*) = 1) which begins at T^ in . In this situation, 



t* = (1— A)t*j+A dt*p(t*)t*. Nonperturbative quark pair creation leading to B — > X c Y(j) 



Figure 5 shows that modest couplings to either hybrids or high mass continua could lead 
to even more substantial duality violations than those associated with the delayed onset of 
the normal cd resonances. 

Although my main focus has been on heavy-to-heavy transitions, the physics issues 
raised here (if not their explicit forms) are also relevant for heavy-to-light transitions. Before 
concluding, let me therefore point out a simple application of the OPE to inclusive heavy-to- 
light transitions where it seems certain to me that they will fail: Cabibbo-forbidden charm 
decays. (Even though such decays might be an unimportant application of the inclusive 
calculations in practice, they provide a valid theoretical testing ground for their accuracy.) 
In particular, consider the c — > diui decays of the D° and D + . They will be dominated by 
the channels D° — > -n^tvn and p~lvg and by D + — > 7r°£z/£, r]£vg, T]'£u£, p°£vt, and uu£u£. Since 
the OPE corrections in the D° and D + are identical, their Cabibbo-forbidden semileptonic 
partial widths and spectral distributions are predicted to be identical. However, simple 
isospin symmetry implies that T(D + — > ir°£ue) = |r(_D° — > n^Iui), so the inclusive Cabibbo- 
forbidden rates can only be equal if T(D + — > i]£i>f) + T(D + —> r^ivy) = T(D + — ► -n^lvi). 
In many models this latter relation would be true if = = m n , since it is rather 
natural for the squares of matrix elements to satisfy its analogue. However, with real phase 




(21) 



where A is the fraction of 



P%n d ue to nonresonant states and pit*) is the appropriate 



may be expected |T{| to persist up to 2 GeV above threshold. 
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space factors, this relation is typically badly broken. Since Cabibbo-forbidden decays, like 
their Cabibbo-allowed counterparts, will receive little excited state compensation given the 
available phase space, I expect this prediction to fail. 

Finally, I note that the duality-violating effects I have highlighted here will have an effect 
on the long-standing B semileptonic branching ratio puzzle [^t|. Since the hadronic mass 
distribution in h — > cud is weighted toward higher masses than the leptonic mass distribution 
in b — > c£i>e, the ratio of these two rates will be changed. 

In summary, I have shown here that hadronic thresholds lead to Aqcd / mQ violations of 
duality in b — > c decays which do not explicitly appear in the operator product expansion. 
Since such violations cannot appear as the b — > c energy release T — > oo, there are "con- 
spiracies" (i.e., sum rules) which relate hadronic thresholds and transition form factors. As 
emphasized by Bigi, Uraltsev, Shifman, Vainshtein, and others JHJ0I, these relations tend to 
compensate the otherwise extremely large Aqcd/^q effects even at small T. In this paper 
I have displayed several models of such hadronic compensation mechanisms which indicate 
that these duality- violating AgcD/ m Q effects could nevertheless be very substantial. While 
the examples I have selected are perhaps pessimistic, they indicate that these effects must be 
better understood before inclusive methods can be applied with confidence to heavy quark 
semileptonic decays. 
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